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The statistical mechanics of a treelike polymer in a confining volume is relevant to the packaging of 
the genome in RNA viruses. Making use of the mapping of the grand partition function of this system 
onto the statistical mechanics of a hard-core gas in two fewer spatial dimensions and of techniques 
developed for the evaluation of the equilibrium properties of a one-dimensional hard rod gas, we 
show how it is possible to determine the density and other key properties of a collection of rooted 
excluded-volume tress confined between two walls, both in the absence and in the presence of a one- 
dimensional external potential. We find, somewhat surprisingly, that in the case of key quantities, 
the statistical mechanics of the excluded volume, randomly branched polymer map exactly into 
corresponding problems for an unrestricted linear polymer. 
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I. INTRODUCTION AND MOTIVATION 

Unlike the structure of the protein envelope of viruses, 
which is well-studied and precisely characterized [TJ [3] 
important aspects of the precise physical organization of 
the packaged genome are as yet undetermined. This is 
true in the case of both DNA [3] and RNA [5] viruses. 
However, it has also been determined that a single strand 
of RNA will organize into a tree-like secondary structure 
[BJE]- Such secondary RNA structure occurs in general, 
and tree-like configurations are known to characterize the 
genomic conformation of certain RNA viruses E]. In 
light of those facts, one can hope to construct a rea- 
sonably accurate theoretical model of the statistical and 
mechanical properties of long segments of single-stranded 
(ss) RNA — and in particular of the genomic matter in ss 
RNA viruses — if one can properly evaluate the statisti- 
cal mechanics of a tree-like polymer in the presence of 
an external potential energy. The potential energy plays 
two roles in the context of the genomic conformation in 
RNA viruses. First, one naturally posits an energetic 
barrier that serves to confine the polymer to a particu- 
lar region in space. Second, given the known interaction 
between RNA and the protein shell (HUE], it is reason- 
able to assume an attractive potential in the vicinity of 
the boundaries of that region. 

The key challenge in this problem is taking into ac- 
count the effects of excluded volume. That excluded vol- 
ume effects are of central importance has been known 
since the late 1940's when Zimm and Stockmayer [T5] 
showed that the radius of gyration R(N) of an ideal 
branching polymer scales with the number of monomers 
N as N 1 / 4 . This scaling relation means that non- 
excluded-volume branching polymers arc highly con- 
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densed; the monomer density grows linearly with dis- 
tance from the center of the polymer coil, a situation 
that cannot be sustained in light of excluded volume con- 
straints. A field theory for excluded-volume interactions 
of branched polymers has been constructed by Luben- 
sky and Isaacson in the form of a 6 — e expansion |13j . 
More recently, Parisi and Sourlas [14] have utilized su- 
persymmetry techniques to argue that the exponents of 
a d-dimensional branched polymer with excluded-volume 
interaction described by this field theory can be obtained 
by a mapping to the Yang-Lee edge singularity [T5J [TB] 
of a d — 2 dimensional Ising model. For d — 3, this leads 
to the scaling relation R(N) cx N 1 / 2 , with a density pro- 
file that now decreases inversely proportional to distance. 
The supersymmetry method is a demanding formalism, 
but Brydges and Imbrie [T71 [TB] showed that it could be 
reformulated as a relation between the conformational 
statistics of a branched polymer with excluded volume 
effects in d dimensions and the statistical mechanics of a 
hard-core liquid in d — 2 dimensions. See also Cardy [TS] 
for a particularly accessible exposition. 

A previous paper by the present authors [20 contains 
an account of the utilization of dimensional reduction 
methods introduced by Brydges and Imbrie to determine 
the conformational properties, in particular the density, 
of rooted trees confined to a finite region and subject to 
an external potential. Although dimensional reduction 
holds in a curved geometry (see appendix [E]) , such as 
the interior of a sphere — the geometry most relevant to 
RNA encapsulation in a viral capsid — it does not appear 
to lead to the kind of fundamental simplification that al- 
lows for the analysis of the effects of interactions between 
the branched polymer and the surrounding walls. Con- 
sequently, our attention was focused on the simpler, but 
still relevant, problem of a branched polymer confined be- 
tween two walls parallel in 3 dimensions. We were able 
to allow for interactions between the polymer and the 
walls, as well as any other one-dimensional external po- 
tential energy This paper provides background to that 
shorter work by filling in important calculational details. 
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It also extends the results reported there. In particu- 
lar, we build on the central-and somewhat surprising — 
outcome that the statistical mechanics of the excluded- 
volumc branched polymer maps onto the statistics of an 
unrestricted chain polymer to investigate both the den- 
sity profile of the branched polymer and the interaction 
between bounding surfaces mediated by it. Given the 
substantial history of research on the statistics of unre- 
stricted chain polymers in the literature(see, for instance, 
[2 1 j ) . the results we present here are not entirely new. 
However, we hope that they will prove stimulating in fur- 
ther investigations of the structure and assembly process 
of spherical viruses. 

An outline of the paper is as follows. In Section [TT] 
we review the consequences of the connection between 
the statistics of an excluded volume randomly branched 
polymer and the statistical mechanics of a gas of hard 
rods in two fewer dimensions. We focus on the simplest 
case of a one-dimensional hard rod gas, namely a gas in 
which the external potential is equal to zero. This is in 
order to develop key formulas and, additionally, to build 
some mathematical intuition with regard to the behavior 
of the more general system in which the external po- 
tential is not constant. Section ITTT1 addresses the means 
of solving for the density and partition function of the 
one-dimensional gas when the external potential varies 
spatially. The method utilized is based on an integral 
equation for the density of the gas derived to Percus [22] . 
We exploit the reformulation of that method by Vander- 
lick, et. al. [53], which we recast into a form suitable 
for a "lattice gas" of one-dimensional hard rods. In Sec- 
tion |IV| the approach is further developed, so that the 
density of the gas, and hence the generating function for 
branched polymer statistics, follows from the solution of 
linear recursion relations. The connection between those 
recursion relations and the Schrodinger equation is de- 
veloped in Section [V] This connection provides the jus- 
tification for the close relationship between the statistics 
of the excluded volume, randomly branched polymer and 
the unrestricted linear polymer chain. 



II. DIMENSIONAL REDUCTION: THE MAP 
FROM THE TREE-LIKE POLYMER IN d 
DIMENSIONS TO THE HARD-ROD GAS IN d 2 
DIMENSIONS 



The principal result that we will utilize relates the 
density, n{x, z) of a hard-core gas in d — 2 dimensions 
to the generating function, z) of rooted, branched 
polymers in d dimensions. Here, we assume translational 
symmetry in all but one dimension of the two systems, 
and x is the one coordinate in the direction along which 
there is any spatial variation. The link between the two 



quantities is expressed in the two relationships [TT] [HE US] 
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E d (x,z) = Y,z N Z N (x) (2.2) 



JV=1 



which tells us that the quantity Z^ contains informa- 
tion concerning the number of configurations of an N- 
monomer branched polymer in d dimensions that is 
rooted at the position x (through ( |2.2[ )) or, alternatively, 
concerning the density at x of an .ZV-particle gas with 
hard core repulsion in d — 2 dimensions (through (2.1 )). 
According to the above equations, the behavior of the 
polymer system with positive fugacity, z, is obtained by 
investigating the mathematical structure of the gas in the 
grand canonical ensemble, but with negative fugacity. 



As an example of the application of Eqs (2.1 1 and (2.2 1 



and to establish some points of reference for the discus- 
sion to follow, we will review the statistical mechanics 
of a gas of rods in a very large one-dimensional interval 
subject to a constant potential energy. 



A. One-dimensional rods in an extended region 
under the influence of a constant external potential 

The constant external potential can be set equal 
to zero. The grand partition function of the one- 
dimensional gas is given by 



Q(L) = £ 
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where a is the hard-core "radius." The factor a in the 
denominator guarantees that each term in the summand 
is dimensionless. In the limit of very large L, the sum will 
be dominated by the N value for which the summand is 
maximum. To locate that term we first exponentiate the 
summand and then re-express the exponent in terms of 
the density p = N/L. Then the summand is of the form 
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where we have made use of Stirling's formula: In TV ss 
Nhi(N/e). The extremum equation that follows from an 
attempt to maximize (2.4 1 with respect to N, and hence 
p, is 
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Then, 



and (2.5 1 becomes 
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The function he h on the left hand side of (2.9) is graphed 
in Fig. [I] The solution to Eq. (2.9) is, formally, h — 




2 h 



FIG. 1: The function, he h , on the left hand side of j^9 l 
graphed as a function of h. 



h(z) 




FIG. 2: The real (soli d cu rve) and imaginary (dashed curves) 
parts of solutions to (2.9 1. 



This tells us that, in the immediate vicinity of z = z c 



h = -l± V2e6z 



(2.12) 



or, from (2.7 1 
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W(z), where W is the Lambert W function [23]. Figure[2] 
shows what the solution looks like, as the fugacity varies 
from a negative to positive values. Note the onset of 
imaginary components to the solution. The departure 
from a purely real solution follows from the minimum in 
the function he h , as displayed in Fig. [l] When z lies 
below this minimum, there is no purely real solution to 

(pi 



The behavior of the solutions in the vicinity of the 
point at which the imaginary part emerges yields infor- 
mation about both the location of this "critical point" 
and about its implications for the statistics of self- 
avoiding rooted trees in three dimensions. The minimum 
in he h is at h — h c — — 1, corresponding to z = z c = — 1/e 
in (2.9 1. Expanding both sides of that equation about 



those critical values we have 
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Solving to lowest order for Sh: 



Sh = ±V2e8z 



(2.10) 



(2.11) 



When z > z c continuity of the solution with z > re- 
quires that we take the upper sign in (2.11 1 ( 2. 13| ) . How- 
ever, when z < z c , the choice of sign is controlled by the 
choice of location with respect to a branch cut in the 
complex z plane, starting at the branch cut at z = z c 
and extending to z = —oo along the negative z axis. 

The fact that the (uniform) density of the one- 
dimensional hard rod gas possess a singularity going as 
(z — Zc)^ 1 / 2 , where z c = —1/e, allows us to extract the 
leading behavior of the large N coefficients in a power 
series expansion of this function of the fugacity |25j . 
Making use of the general result that if a function f(z) 
can be written as a power series about z = 0, so that 
/0) = Y,n=o F nz N , then 
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where the contour over which the integral in (2.141 is 
performed encircles the origin in the complex z plane 
and does not enclose any singularities in the function 
f(z). Expanding the contour so that it impinges on the 
branch point and, ultimately surrounds the branch cut 
(see Fig. [3]), we end up with the following result for the 



4 



coefficient of z N in the power series expansion of p(z): 
1 r°° ( 
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(2.15) 



According to (2.1l and (2.2), we see that the num 



ber of conformations of three dimensional trees with AT 
monomers should grow geometrically as l/(irz c ) N , with 
the additional power law modification N~ z / 2 . 

Im (z) 



Re (z) 



FIG. 3: The original contour in ( |2.14[ ) (small circle surround- 
ing the origin) and the distorted contour that bounds the 
branch cut originating at the branch point, z c . 



Note that in the above discussion, we have neglected 
the possibility of any other non-analyticity in p{z). As 
one can readily verify, the leading contributions to the 
integral in (2.14) will, for large N, be controlled by the 



singularities in p(z) that lie closest to the origin. The pos- 
sible existence of other poles, branch points or essential 
singularities in p(z), all of which will occur at \z\ > \z c \, 
is irrelevant to the results obtained here. 

The fact that the statistics of the hard-core gas — and 
by extension randomly configured trees — are controlled 
by non-analyticity in the grand canonical ensemble that 
is closest to the origin in the complex z plane can be ex- 
ploited to infer the emergence of a bound state generated 
by the presence of an attractive external potential. The 
argument is as follows jTS] . 

We start with the density at a particular point of the 
gas of hard rods. In the grand canonical ensemble, this 
density equal to the sum over configurations in which 
the a rod is at that point, divided by the sum over all 
allowable configurations. If the system is confined to the 
one dimensional region between x = and x = L, then 
this leads to 



p(x,z) 



Q(x, z)zQ(L 



a) 



Q(L,z) 



(2.16) 



with Q the grand partition function as before. If, now, 
there is a non-zero, delta function potential energy, V 
at the position x±, then the grand partition function 
contains an additional contribution from the associated 
Boltzmann factor weighting configurations in which a gas 
particle occupies that particular location. If we write 



~ fJV = 1 + A 



(2.17) 



where an attractive potential (negative V) generates a 
A > 0, then the new partition function is 

Q'(L,z) = Q(L,z) + Q(x u z)\zQ(L-xi,z) 

= Q(L,z)(l + \p(xi,z)) (2.18) 

The new density at x\ is, then, given by 

Q(xi,z)zQ(L — xi,z) 



p'(xx,z) 



Q'(L,z) 
Q(x 1 ,z)zQ(L - x x ,z) 
Q(L,z)(l + Xp( Xl ,z)) 



(2.19) 



In the case of a very long interval and a large value of 
X\, the density p{x\, z) will be independent of x\, and is 
given by (2.5 1-( [2~9| . Figure [1] illustrates the behavior of 
the density as a function of z for z > z c — — 1/e. The 
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FIG. 4: The function p(z), as given by (2.5 1-12.9 1, for z > z c , 



the location of which is indicated by the vertical dashed line 



fact that p(z) is negative for 

— oo as z — » 



< and that it goes to 



(see also (2.13 1 ) ensures that, for any 



positive value of A, there will be a zero in Q'(L.z) and 
hence a pole in p'{x\, z) for a (negative) value of z closer 
to the origin than z c . If we associate this singularity in 
the density with a bound state, we are led to conclude 
that an arbitrarily weak one-dimensional potential well in 
the interior of an infinitely extended system will "bind" 
a three dimensional random tree. Further analysis of 
this situation requires a more searching exploration of 
the spatial structure of this putative bound state. 
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III. ONE DIMENSIONAL HARD CORE GAS IN 
THE PRESENCE OF AN EXTERNAL 
POTENTIAL. 

The problem of a one dimensional gas of hard core 
rods was considered by Percus |22j . who derived an in- 
tegral equation for the equilibrium density of that sys- 
tem, from which quantity the partition function and all 
interesting equilibrium correlation functions can be ob- 
tained. The integral equation, which has been reduced 
to a very useful and tractable form by Vanderlick, et. 
al. |23) , serves as a starting point for the exploration of 
the equilibrium statistics of randomly branched polymers 
in a one-dimensional environment , given the connections 
established by Brydges and Imbrie [TTJ [TS]. The two 
equations leading to the calculation of the density in this 
one-dimensional system are [22j 123] 



h(x) 



h(x) 



p(x) 



i-r x +a P (t)dt 



e 0(p-u(x)) exp 



h(t)dt 



(3.1) 



(3.2) 



where (3.1) can be taken as a definition of the function 
h(x). 

Note the strong similarity between (3.2 1 and (2.9 1. In 
fact, if we set u(x) = 0, replace e 13 ^ by z, t ake a = 1, 
and assume an ^-independent h(x), then (3.2 1 reproduces 



(2.9), while, if we also assume an ^-independent density, 
p, and set a = 1, then (3.1l reduces to (2.6 1. 



For our purposes, it proves more useful to focus our in- 
vestigations on a discrete version of the one-dimensional 
gas of rods — a one-dimensional, hard core lattice gas. As 
it turns out, the equations leading to a solution for the 
density of this system have already been worked out in 
a different context |26) . To maintain a self-contained ex- 
position, the derivation of the equations governing the 
density in this discrete system are also presented below. 
This derivation closely parallels the arguments of Percus 

We assume rods that can sit on specific points on a 
line. As shown in Fig. [5] the possible locations for the 



the gas lie between the sites at the points indicated by 
the x 's. These boundaries are utilized in the definitions 
of "partial" grand partition functions. For instance the 
grand partition function S(— oo, n) corresponds to a sys- 
tem in which the rods can occupy the region to the left 
of the x at n in Fig. [5] The full grand partition function 
of the infinite system is S(— oo,oo) = 3. Now, assume a 
chemical potential p and site-dependent local potentials 
u n . We locate a rod by indicating the site at which its 
far right portion sits. For instance, Fig. [6] shows a rod 
with a length equal to two that is located at the site n. 
Given this, it is straightforward to show that the density 



n-1 n n+1 



qj_l 



FIG. 6: A rod with a lenth a — 2 located at the site n. 
of rods at the site n is given by 

5(-OQ,Tl- (7)5(^,00)6^-"") 
Pn = ^ (3-3) 

The next step is to establish a relation for the product 
of two partial grand partition functions. In particular, 
we are interested in S(— oo, n)S(n, oo). This product is 
almost the entire partition function. In fact, it only omits 
configurations in which a rod overlaps the boundary at n. 
In the case of the rod pictured in Fig. [6] there is precisely 
one such configuration: the one in which the right hand 
side of the rod lies on the site n+1. In general there 
are a — 1 such configurations. Let us look at the one 
corresponding to the location of the rod with length a 
lying on the site The contribution to S of that 

particular configuration is 



S(— oo, n + 1 — cr)H(n + 1, oo)e 



= So 
^Pn 



(3.4) 



Making the appropriate corrections for all missing con- 
tributions to S, we see that 



n-1 n n+1 



n-1 n n+1 



FIG. 5: The setup for the one dimensional hard rod gas. 



rods are at sites labeled by the integer n. The actual 
locations of the sites are indicated by the vertical lines in 
the figure. The boundaries between different regions of 



E(-oo, n)^(n, oo) = H x ( 1 p n+j 



(3.5) 



Now, we construct a recursion relation for the partial 
grand partition function. By inspection, one sees that 
the difference between the partial function S(— oo, n) and 
the function S(— oo, n — 1) lies in configurations in which 
there is a rod at n. In other words 

S(-oo,n) = S(-oo,n-l) + e /3 ^- u " ) S(-oo,n-cr) 



= S(— oo, n — 1) 
= S(— oo, n — 1) 



S(n, oo) 
/9„S(-oo,n) 

1 ~ Ei=l Pn+j 



(3.6) 



6 



The second line of (3.6 1 follows from (3.3) and the last quantity h n , given as 



line from (3.5 1. Thus 



H(— oo,n — 1) = H(— oo,n) 



Pn 



1 - Z)J=1 Pn+j 



a(-oo,n)(i - /i„) 



(3.7) 



1 _ X^=i Pn+j 



(3.8) 



where the last line of (3.7l serves as a definition of the (3.3), we note that, given (3.7l 



Returning to the expression on the right hand side of 



_ Un) S(-oo, n)(l - h n )(l - feg-i) ■■■(!- fe»_ g+ i)^(n, oo) 



j — 1 m— 1 



(3.9) 



We then divide both sides of (3.9) by the term 1 



ST=i Pn+j- Then, making use of the definition of h). 



implicit in the last line of (3.7 1, and singling out the first 



term in the product on the right hand side of (3.9 1 we 
have 



Another version of the set of equations that is a bit 
simpler to iterate replaces the variables h n by 



9n = 1 - K 



Then, 



(3.14) 



h n = e ^-«»)(l - h n ) H (1 - K- m+l ) (3.10) 



Solving for h n , we are left with 



h n = 



nl=2( 1 _ hn-m+l) 



l + e /3(^)n^ 2 (l-^ m+1 ) 



If we define 



H„ — Y\_ (1 hfi-m+l) 

m=2 



then the equation (3.111 becomes 



1 + H n eP(»- u *) 



(3.11) 



(3.12) 



(3.13) 



The quantity H n can be thought of in terms of a local 
modification of the fugacity arising from the excluded 
volume constraint. Note that the quantity h n is deter- 
mined by h n 's to the left of it. This means that the 



system of equations (3.11) can be solved by iteration to 



the right. Furthermore, given the explicit definition of h n 
in equation (3.8 1, we see that the density, p n is expressed 
in terms of h n and p m 's with m > n, so having solved for 
h n , we obtain the p„'s by iterating to the left. 

To highlight the precise points of reference between 
the discrete equations above and the integral equations 
developed by Percus [22] an d refined by Vanderlick, et. 
al. |23j. we note that (3.8 1 corresponds to (3.1l while 



(3.111 is the analogue in the discrete system of (3.2 1 



H n — 9n-m+l 

and the equation for g n becomes 

1 



9n = 



1 + Hr.e^-^) 



with 



_ , Pn 

In — 1 ^<r-l 

1 _ 2^7 = 1 Pn+j 



(3.15) 



(3.16) 



(3.17) 



The equations above constitute discrete version of the 
method of Vanderlick, et. al. |2"3"1. 



A. The case a = 1 

It is worthwhile to ask what happens when a = 1. In 
this case, the rods are "point particles," and excluded 
volume does not play a role. The solution of the equa- 
tions is greatly simplified. Given (3.9 1 and (3.8 1 with 



(7=1, one arrives at the following end-result. 



Pn 



e /30-«n) 
1 _|_ e P(n-u„) 



(3.18) 



The lattice gas in this case is, in fact, well-described by an 
appropriate version of the one-dimensional Ising model. 
This system maps onto trees in which the "branches" 
have zero extension, and in which, furthermore, the ra- 
dius of the hard core at the vertices between branches is 
vanishingly small. 
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B. The lattice gas with a — 2 

The simplest non-trivial version of this lattice gas as it 
applies to the conformational statistics or random trees 
has a = 2. Then (3.11 I and (3.8 1 become 



h n = 



1 



S^-"")(l-ftn-l) 

fe^ifl-^-!) 

Pn 
- Pn+1 



(3.19) 
(3.20) 



For the time being, we will focus on a uniform gas 
confined to a finite region. This means that we are going 
to set u n = 0, supplemented by boundary conditions to 
be developed below. Replacing by the fugacity, z, 
the equation for h n becomes 



h n = 



z(l - h n -i) 
1 + z(l - /i„_i) 



(3.21) 



In terms of g n , as defined in (3.141, (3.211 becomes 



1 



l + zg n -i 



(3.22) 



This equation leads to a solution for g n in terms of g\ in 
terms of continued fractions: 



1 



1 



(3.23) 



We can rewrite the solution to the equation for g n in 
terms of an iterated matrix equation |27j . Given the fact 
that the continued fraction is repeated, the iteration of 
the matrix equation is relatively straightforward. Sup- 
pose we have the solution for g m in terms of g m -i as 
follows: 



Then, given 



we have 



9n 



9m-l 



At + Big 



m—l 



Ci + Dig m _i 



1 + zg m -i-i 



(3.24) 



(3.25) 



9i = 



Ai+i + Bi +1 g m -i-i 
Ci+i + -Df+lPm— /— 1 

Ai + Bi t—± 

1 « l + ZO„_,_ 1 



(3.26) 



Rationalizing the last line of (3.261, we end up with the 
following recursion relations 



Equations (3.271 and (3.28) can be written as follows 

(3.31) 



Ai+i 
Bi+i 



1 1 

z 



Ai 
B, 



with a similar equations for C and D. Given that we 
have the initial conditions Aq = 0, Bo = 1, Co = 1 and 
Dq = 0, we have 

i 

(3.32) 



.4, 
B t 



1 1 

z 



Our task now is to take the I power of the matrix in 



(3.311 



C. A digression: the uniform case 

Before performing the requisite calculations, we will 
consider the solution (3.221 under the assumption that 



the g's are independent of location. The equation that 
results from that model is 



-1 = 



The solution is 



(3.33) 
(3.34) 



There is a branch point at z = —1/4. Figure [7] shows 
what the real and imaginary parts of g look like as a 
function of z. Close to the critical point, the general 



h(z) =l-g(z) 




FIG. 7: The real (solid curve) and imaginary (dashed curves) 
parts of h(z) = 1 — g(z), where g(z) is as given by p34| . 

form of the plot is qualitatively identical to the result for 
h(z) displayed in Fig. [2] 

We can reconstruct the (constant) density, p, from 
(3.201, with h = 1 — g, where g is given by (3.34). From 



(3.201 with p, h and g constant, we have 
h 

1 + h 



Ai+i 


= A+B t 


(3.27) 


1 - 


g 


Bi+i 


= zAi 


(3.28) 


2 - 


g 


Ci+i 


= Ci + Di 


(3.29) 


2z 


f l 


A+i 


= zBi 


(3.30) 


1 + 


4z 



4z - VI + 4z 



(3.35) 
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If z = -1/4 + 6, then 



P 



1/2 + 25-2^/5 



4(5- 
1 

~4v^ 



2VS 



(3.36) 



The solution (|3.34|) suggests a reparameterization that 

(3.37) 



is useful in the vicinity of the "critical point." 

1 



4 cosri k 



Then, (3.341 becomes 

g = 2coshke- k (3.38) 
When z passes through -1/4, (3.371 requires that we re- 



place the real k by an imaginary quantity. 



n 



9! 

4 ■ 

3 - 
2-^ 



50 



100 



150 



200 



n 



FIG. 8: The function g n when z is slightly to the right of the 
critical point at -1/4. The dashed line at g n = 2 indicates the 
value to which g n tends when z = —1/4 exactly. 



D. Completion of the calculation 

The results we need follow from the eigenvectors and 
eigenvalues of the matrix ( ^ ^ V Solving the rele- 
vant equations and making use of the reparameterization 
(3.371, we find that (3.321 reduces to 



. \ / sinh kn \ 

A n \ _ / (2coshfc)"- 1 sinhfc I 

I I sinh(n— l)fc I 

" / \ ~~ (2coshfc)™ sinhfc / 

With a similar set of steps, we end up with 

sinh(n+l)fc 



D. n 



sinh k{2 cosh fc) 11 
sinh kn 



(3.39) 



(3.40) 



\ sinh k(2 cosh fc)' 1 + 1 / 

We assume that the initial g is equal to zero. This 
is consistent with assumptions in the case of the ordi- 
nary one-dimensional excluded- volume gas. We will as- 
sign this g the index 0. Then, 

gn(z) = -^r 

2 sinh nk cosh k 



sinh(n + l)fc 



= 2cosh 2 fc-sinh2A:coth(n-(- l)k (3.41) 

If k is real, the n —> oo limit of the above expression is 
consistent with (3.35). However, if k is imaginary, then 



the second term on the last line of (3.41 1 behaves like the 
cotangent function, and is periodic with a set of poles. 
See Fig. [HJ The graph in this figure was g enerat ed by 
numerical solution of the recursion relation (3.221. The 



key property of this plotted solution for g n , and in the 
analytical form in (3.41 1, is the appearance of poles. In 



fact, as will be demonstrated below, the passage of a pole 
in g n through the far boundary of the interval of interest 
is directly associated with the appearance of a singularity 
in the dependence of the density on the fugacity, z. A 
fuller discussion of this point will follow the development 
of an alternative approach to the solution of the recursion 
relation for g n (z). 



IV. MODIFIED APPROACH TO THE 
SOLUTION OF THE RECURSION RELATIONS 

We now describe an alternate route to the analysis of 
this system. Recall Eq. (3.221. Suppose we write 



1 + zg n (z) 



(4.1) 



Then, (3.221 is transformed to the following recursion 
relation 



= 1 



which is manipulated to 

4>n+l(z) = 1p„(z) + Zlj) n -i(z) 

Now, let 



(4.2) 



(4.3) 



Mz) - (-z) n/2 0„(z) 



1 

2 cosh k 



4>n(z) 



(4.4) 



where in the last line of ( 4.4 1 use has been made of (3.37 1 



Substituting from (4.4 1 into (4.3 1, we end up with the 
recursion relation 

(j) n+ i(z) - 2coshfc0„(z) + (j) n -i(z) = (4.5) 

The solution to this equation is 

(j) n {z) = Ae kn + Be- kn (4.6) 

In the case that k is imaginary, corresponding to z < 
— 1/4, the solutions are complex exponential — or, replac- 
ing k by ik, 



4> n (z) = Acqs kn + B sin kn 



(4.7) 
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A. Boundary conditions 

We require that the density be zero at the two ends of 
the interval. We set the leftmost end at n = and the 
rightmost one at n = L + l. Making use of (3.8) with h n 
expressed in terms of g n , we have 



Pn = (1 - <7n)(l ~~ Pn+l) 



(4.8) 



We ensure that po = by requiring go = 1. As for 
Pl+i, we simply set it equal to zero. The pole in the 
solution for the density will be due to a pole in g^ , which 
will, according to (4.1 I follow from iPl(z) passing through 
zero. The boundary condition at n — will be obtained 
if we set 



ipo(z) = ipi(z) = 1 



(4.9) 



Making use of the above equation and Eq. (4.1 1 we see 
go = The pole at n = L follows from the requirement 
that 



tf> L (z) = <f> L {z) = 



(4.10) 



Given that k will be small in most of the cases of in- 
terest to us, our solution in the uniform case will be gen- 
erated from 



4> n = sin(Z7r(n + l)/(L+l)) 



(4.11) 



with I an integer. Here is how the argument goes. We 
start with (j> n of the form sin(kn+8) . Then, given that the 
ratio tpi/tpo is equal to one, making use of the relationship 
(4.4) between tp n and <j> n , we have 



sin(fc + 9) 
2 cos k sin 9 



1 



(4.12) 



Anticipating that both k and 9 will be small, we expand 
the sine and cosine functions in terms of their arguments, 
and we end up with 



k + t 



= 1 



(4.13) 



The solution to this equation is 9 = k. Requiring that 
kL + 9 = In, we end up with 



lw/(L + l) 



(4.14) 



Substituting the results into the argument of the sine 
function, we end up with (4.11). 



B. Reconstruction of the density: general formulas 

The next step will be to reconstruct the density from 
the function g n (z). We can do this formally, by noting 
that the general form of the equation (3.201 is 



Pn+i 



PnPn 



(4.15) 



We can solve this equation by analogy with the solution 
of a linear differential equation. We start by introducing 
a new variable 



Rn — Pn ] i Prn 



(4.16) 



Making use of this new variable, (4.15 1 is recast into the 
recursion relation 

n 

Rn+l = Rn + OL n J~J Z?" 1 
m— mo 

= Rn+ln (4.17) 

The solution of this very simple recursion relation is 

n-l 
l=l 



In the case of the product in (4.16), the understanding 



is that if n — 1 < tuq, then the product is re place d by 
rim=n-i A" and simil arly for the product in (4.171. As 
for the sum in (J4.18 1, if n — 1 < Iq then the s um i s 



replaced by — Y^i= n H- Finally, if n — 1 = mo in (4.16 1, 
then the product is replaced by /3„_i =mo , and similarly 
for the product in (4.171. Note that all considerations 



are simplified if we take mo = Iq = 0. 

Ultimately, we recover the density p n through the in- 
verse of (4.161, with one proviso. Taking m = l a = 0, 



we note that a full solution to (4.15) is of the form 



f n-l \ n-l 
\ 1=0 / m=0 



(4.19) 



with C an arbitrary constant. In the case of interest, the 
equation corresponding to (4.151, as given by ( 3.20[ ), is 



so here, 



Pn+l = 1 - h n X p n 



a n = 1 

Pn ~ h"n 



(4.20) 



(4.21) 
(4.22) 



C. Specific relations for the reconstruction of the 
density 



Given the relationships (4.1) and (4.4 1, we can write 



Solving for g n (z), we have 

o (z) - -4cosh 2 k ( 1 

g n [z) 4cosn ^ 2cQshk 



- 1 



(4.24) 
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Now, we use (4.111, with k general and assume trigono- 
metric functions (k — > ik). Then, 



(j) n {z) = sin k{n + 1) 



(4.25) 



Where, recall, z = — 1/4 cos 2 fc. Substituting this into 



(4.241, we have 

hn(z) = 1 - g n (z) 

= 1-4 cos 2 k + 2 cos k 

sin k{n — 1) 
sin kin + 1) 



sin k(n + 2) 
sin k(n + 1) 



(4.26) 



As our next step, we will implement the recursion re- 
lation repeatedly. Starting with a particular value of n, 
which we will set equal to L, and assuming pl+i = 0, we 
have 



PL 



sin k(L — 1) 



sinfc(L+ 1) 
sin k(L — 1) 
smk(L + 1) 



(PL + l - 1) 



(4.27) 



If we apply this recursion relation down to pl-3, we end 
up with 



PL-3 



sinfc(i-4) s'm k(L - 4) s'mk(L - 3) 
~sinfc(i - 2) ~~ sinfc(L- 2)sinA;(L- 1) ~ 
sin k(L - 4) sin k(L - 3) sin k(L - 4) sin k{L - 3) 



sin fc(L — 1) sin k(L) 



sin k(L) sin fc(L + 1) 

(4.28) 



There will clearly be a pole in the right hand side of 
(4.281 when sinfc(L + l) = 0. However, the denominators 
in other terms that expression also pass through zero, at 
different values of k. As we will see in our discussion of 
the more general case, those denominators do not give 
rise to additional poles. 



D. The general case 

We now turn to the way in which the reconstruction 
of the density works in the most general case. That is, 
we look into what happens when the fugacity varies from 
position to position. In that case, the recursion relation 
can be written in the form 



1 



1 + z n 



If we write 



ZjiQfl 



lb. 



n+1 



(4.29) 



(4.30) 



we obtain the following relationship: 

lpn+2 = i>n+l + Z n+ i1pn 



Then, 



h n = 1 

= Z n -1 



On 

Ipn 



(4.32) 



We can reconstruct the density in the same way as we 
did in the case of no potential. Carrying out the same 
procedure as above, we find 



PL-3 



zl-4- 



L-5 



+ZL-4Z L - 3 Z L - 2 



Z L _4Z L _3- 



L-4 



-Z L -4Z L - 3 Z L - 2 ZL-1- 



ib L _ 3 ib L _ 2 

1pL-2lpL-l 

ib L - 5 ib L -4 



(4.33) 



Once again, the pole in the density results from a zero 
in tpf To see that no other zeros lead to a pole, we 
consider the case of ipL-2- The terms in (4.331 in which 
that function appears combine as follows. 



1 



1 



ZL-2 



-ZL-4ZL-3lpL-5lpL-4 
-Z L -4ZL-3lpL-5^L-4 



V>i-2 \1pL-3 i>L-\ 

1 ( 1pL-l - Z L -2lpL-3 



1 lpL-2 
^L-2 l\)L-\i>L-3 



(4.34) 



where the last line of (4.341 follows from (4.311 with a 



suitable adjustment of n. 

We can alter the equations for ip m and p m by intro- 
ducing a modification of the function 4> n . Our new and 



generalized version is via the following alteration of (4.4 1 



Then, (4.311 becomes 

{-Z n+2 ) 1/2 4'n+2 = 0n+l 



(4.35) 



(4.36) 



Furthermore, the last term in ( |4.34 1, the term in which 
there is a pole resulting from a zero of ipL — an d hence 
4>l — reduces to 



{-Z L -i) 



1/2 



>L-5<PL-4 



(4.37) 



RELATION TO THE SCHRODINGER 
EQUATION 



The recursion relation (4.5), and particularly the solu 



tions (4.6 1 and (4.7 1, strongly suggest a relationship be- 



(4.31) tween the equations that we solve for the quantity g n {z) 
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and thence the density p n (z) and the Schrodinger equa- 
tion for a free particle. In fact, when the variable k is 
small and imaginary, Eq. (14. 5| becomes 



<t>n+i{z) - 20„(z)(l - k 2 /2) + 4>n-i{z) = (5.1) 

which reduces, under the assumption of slowly- varying 
4> n (z), to 



d 2 (t>n{z) 

dri 2 



k 2 (j) n (z) 



(5.2) 



Now, consider (4.361. We write 



1 



4 cos 2 k 



-/3u„ 



Expanding in k and assuming that both 4> n (z) and u n 
are slowly varying functions of n, this recursion relation 
becomes 



d 2 (j)n{z) 

dn 2 



+ (e*- - 1) <f> n (z) = k 2 K(z) (5.4) 



If (3u n is small, the second term on the left hand side of 



(5.4) is just the standard potential energy contribution 



to the Schrodinger equation — with the multiplicative fac- 
tor [3. Under the conditions described above, the recur- 
sion relations and, more particularly, the density through 
(4.371, are obtained via the solution to the Schrodinger 



equation. In fact, again, if (3u n is small, the residue of 
the pole in the density as a function of z is given by 



</)„(z) 2 Res(l/0 L _i(z)^( 2 )) 



(5.5) 



In fact, the residue function can be directly related to 
the normalization of the solutions to the equation (5.4 1. 



This fact can be established by appealing to a standard 
result for the normalization of the solutions to the one- 
dimensional Schrodinger equation |28j . This argument 
is contained in Appendix [X] For an outline of a demon- 
stration based on the discrete equations, see Appendix 

m 



VI. ALTERNATE DERIVATION OF THE 
SCHRODINGER EQUATION FORMALISM 
FROM A GRADIENT EXPANSION 

The fact that we are led by the developments in Sec- 
tions II VI and [V] from a discrete version of the one- 
dimensional lattice gas problem as formulated by Percus 
and Vanderlick to the continuous Schrodinger equation 
suggests the possibility of a direct route from the con- 
stitutive equations of the continuous hard rod gas to the 
Schrodinger equation approach. In fact, guided by what 
has been described previously in this article, one can out- 
line just such a development, based on a gradient expan- 
sion of (3.1 1 and (3.2 1. We start with an analysis of the 



e !3(p-u(x)) exp 



second of those two equations 
h(x) = 

= z(x) exp 

where 



h{t)dt 



h(t)dt 



z(x) = ze~^ 



(6.1) 



(6.2) 



When the system is uniform, the singularity occurs on 
the vicinity of z(x) = zq = — l/ere, h(x) = ho = —X/a. 
We rewrite the (6.1 1 as follows 



(5.3) crh(x) 



:> —crh(x) 



a(-l/(i + 5h{x))e a[ - 1/,J+Sh{x)) 



1 



= -e -i (l - aSh(x))(l + aSh(x) + -a z 8h{xy + ■■■) 



-e- 1 (l + ^a 2 Sh(x) 2 + ---) 



a(-l/ae + Sz)e-^ u ^ 



exp 



h(t)dt - ah(x) 



-e-\l ~ ae8z)e-^e' j2dSh ^' dx ' 2+ - 
-e _1 (l - aeSz)(l - (3u(x) + ■■■) 
,„ 1 9 dSh(x) . 



(6.3) 



Equating the fifth and eighth line of (6.3 1, and expanding 



to second order in a and first order in f3u(x), we end up 
with the following equation 



l- 1 -* 2 5h(x) 2 = l-ae5z + 1 -a 2d6Hx) 



dx 



(3u(x) (6.4) 



dSh(x) 



dx 



- 6h{x) 2 + -^Pu(x) 



2e , 



If we set 



Sh(x) = -W(x) 



then ( |6.5[ ) becomes 
dW(x) 1 



dx 



2 a z 



4e 



(6.5) 



(6.6) 



(6.7) 



Let 



W(x) 



ip{x 
d 

dx 



Then, 



In [V>(x) 2 ] 

d 2 i>(x)/dx 2 



(6.8) 



(6.9) 



12 



and (6.7 1 becomes 

d 2 4>{x) 



dx 2 



2 



2e 



Sz i>{x) (6.10) 



We next turn to the other constitutive equation, (|3.1 I 
To analyze this relation, we multiply both sides 
and perform the same kind of expansion as we did on 
(6.1). To facilitate this expansion, we replace p(x) by 
—p(x), anticipating that the density becomes negative 
in the regime of interest. In fact, it will become large 
and negative, as we are interested in the behavior in the 
vicinity of a pole. Performing the gradient expansion 
utilized above, (3.1l becomes 



a(-l/a + 5h(x))(l + ap(x) + -a 2 



dp{x) 
dx 



+ •••) 



-(1 - a5h{x)){\ + ap{x) + I^^M 

2 dx 



+ ■■•) 



-(1 + <jp{x) — a 2 p(x) Sh(x) 

-aSh(x) H ) 

-ap(x) 



1 2 dp(x) 
2 a dx 



(6.11) 



The neglected terms in (6.111 are higher order in a. 
Equating the last line in (6.111 to the next-to-last line, 
we end up with the equation 

1 



^dp^x) _ a 2,, §h ,, = adh ,, 

2 dx 



1 



Making use of (6.6), this equation becomes 
dp(x) 



dx 



W(x)p(x) = —(aSh(x) - 1) 



(6.12) 



(6.13) 



Note that the solution to the homogeneous version of 
([6~ri is 



p{x) 



exp 



W(x')dx' 



exp [ln(^(x) 2 - iP(x ) 2 )] 
A^{x) 2 



(6.14) 



where we have made use of (6.8 1 



To further analyze (6.131, we replace 5h( x) o n the right 
hand side by tp'(x) /ip\x) , as mandated by (6.6 1 and (6.8 1. 
If we further set 



p(x) — tp(x) 2 r(x) 



Eq. (6.13) becomes 



tp(x) 



2 dr(x) 



dx 



2 / dip(x)/dx 
a 2 \ 4>{x) 



(6.15) 



(6.16) 



In keeping with our gradient expansion approach, we fo- 
cus our attention on second term in parentheses on the 



right hand side of (6.161. The equation that results from 



ignoring the first term is 



x dx' 

o 



3y a and we obtain the following expression for p(x) 



P( X ) = 2^ X ) 2 



We now note that the function 



x dx' 



V>(x') 2 



(x) = <ip(x) , _ 



x dx' 

o W) 



(6.17) 



(6.18) 



(6.19) 



is also a solution to Eq. (6.101, and that it has a Wron 



skian of one with the function ip( x )j m that 4>'(x)rp(x) — 
ip'(x)4>(x) = 1 [29 . Both of these properties can be ver- 
ified by direct substitution. Note that they are indepen- 
dent of the lower bound of integration, Xq. A general 



version of (6.19) is 



p(x) = ip(x)(<j)(x) + Aip(x)) 

a 



(6.20) 



We will assume that ij}(x) satisfies the boundary condi- 
tion ijj(0) — 0, which yields a p(x) that is also zero at 
that interval boundary. To insure that p(L) — 0, we set 
the constant A equal to —<fi(L)/ip(L). There is then a 
pole in the density when ip(L) = 0. To determine the 
residue at that pole, we note the following 



a 



2 nx) - 



2 2 mm 



-MX) 



ip(L) ip'{L) 



(6.21) 



The last line of (6.21 1 follows from the Wronskian rela- 
tion between 4>(x) and ip( x )- The pole term on the right 
hand side is the one going as (if)(L)4>' (L))^ 1 . Follow- 
ing the reasoning in Appendix |A| (see especially (A4|) 
we find that this generates the factor —a/(2e)l/((Sz — 
Szo) Jq tp(x) 2 dx), where 5zq is the value (actually, one of 
the values) of 6z at which i()(L) = 0. Thus, the residues 
at the poles incorporate the normalization of the eigen- 
functions there, as was the case in the discrete version of 
the model. 



VII. APPLICATION OF THE SCHRODINGER 
EQUATION FOMALISM: THE CASE OF A 
DELTA-FUNCTION POTENTIAL 

In a case of particular interest to us the potential en- 
ergy, u n , is non-zero at one site, uq. Here, it is not correct 
to treat the potential energy as either always small or 
slowly-varying. However, as we will see, the connection 
between solutions to the Schrodinger equation and the 
density as constructed from discrete recursion relation 
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still holds. The recursion relation to which this potential 
energy leads is 



9n + l -. _o 

l + ze P u g no 
1 



1 + z '9n 



(7.1) 



Rewriting the g's in terms of <j)'s, we end up with the 
relation 

4>n +2 
= 2 COS k(f>n +X ~ 4>n 

z'\ 2 COS k(f> ng - (f>n + l 



+<t>n 1 - - 



z J 2cosfc(l - z' /z)4> no + 4> na+1 z' /z 



^n +2 



1 - 



z J 2cosfc(l - z'/z)(f> no + 4> no+ iz'/z 

(7.2) 



The quantity <$L +2 1S the value of 4> nQ +2 in the absence 
of the delta function potential. If that potential is small, 
then (1 — z' I z) will be small, and the right hand side of 
(7.2) becomes 



4>n +2 + <t>n - l) 4>no-\/<t>n +l 



(7.3) 



Now, we once again take the ratio with which we 
started. This leads us to the equation 



t>n +2 _ Vn a +2 



t>n a + l <Pn + 



2_ + fZ_ x \Mn^ 

i \zf > <j>i 0+1 



4) 



We will consider two cases: extended states and the pos- 
sibility of a bound state. The latter solution exists if the 
potential is attractive and exceeds a threshold value. 



1. Extended state solutions 

The assumption that we now work with is that the 
solutions to the equations have the following form 



asin(fc(n + 1)) 



n < no + 1 



6sin(fc(n + 1) + 9) n > n + 1 



(7.5) 



Note that there are two possible forms for of <\> n at n — 
no+1. That is, asin(fc(ri +2)) = bs'm(k(n +2)+6). The 
quantity 0° o+2 is, a ccord ing to (7.5 1, equal to asin(fc(n+ 
3)). Given all this, (7.4 1 becomes 



sin(/c(n + 3) + 9) 
sm(k(n + 2) + 6) 
sin(fc(no + 3)) 
sin(/c(n + 2)) 



- 1 



i(k(no + 1)) sin(fcno) 
sin(/c(n + 2)) 2 



(7.6) 

We now consider the most general version of (7.6) as a 



condition on the phase shift 9. That is, we look at the 
equation 



sin(fc(n + 3) + 6) 
sin(fc(n + 2) + 9) 



A 



(7.7) 



Expanding the sine functions and solving for 9, we end 
up with the following result: 



sin(fc(n + 3))-Asin(/c(n + 2)) 

tanfr=— j—, —rz -7—. ~rr" (7. 

Acos(fc(n + 2)) - cos(fc(n + 3)) 



Inserting the right hand side of (7.6 1 into (7.8) as a sub- 



stitute for the quantity A, we end up with the result for 
tan 9: 



tan 9 = 



(z/ z 1 — 1) sin(fc(no + 1)) sin(fcno) 



sin k + (z/z' — 1) sin(/c(n + 1)) sin(fcn ) cos(fc(n + 2))/ sm(k(n + 2)) 

I 



(7.9) 



We are now in a position to work out the allowed values 
of the quantity k. Given the boundary condition (4.101, 
the requirement on k is 



sm(k(L + 1) + 6(h)) =0 



(7.10) 



origin than —1/4. The solution to the recursion relations 
for (f> n (z) will then be 



, asinh(fe(n + 1)) n < uq + 1 
l ^ Z > = ' ber kn n > n + 1 



(7.11) 



2. The bound state 

Here, we assume hyperbolic functions, corresponding 
to a form for k that places the fugacity, z, closer to the 



This solution is appropriate to a system in which the 
length of the region to which the rods are confined is 
arbitrarily great. Applying (7.4 1 to this conjectured so- 
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lution, we end up with the relationship 

sinh(fc(n + 3)) 
_ sinh(fc(n + 2)) 

/ z \ sinh(fc(n + 1)) sinh(feroo) 
+ V?~ ) sin 2 (fc(n + 2)) 

which can be manipulated to 

sinh fcsinh(fc(no + 2)) 
= (l - J) sinh(fcn ) sinh(fc(n + l))e~ k ^ 



(7.12) 



(7.13) 



Limiting cases 



Two limits are of interest to us. First, if hq 3> 1, then 
(7.131 reduces to 



2 sinh fee 



(7.14) 



Given I — z/z' = 1 — e' 3 ", we note that (7.141 has a 
solution of the type desired only if u < 0. If we expand 



the left hand side of (7.141 in u and, assuming negative 



u, replace the left hand side of (7.14) by — \(3u\, we obtain 
the following equation for k 



„4fc 



'2k 



\f3u\ = 



The solution to this equation is 
1 



k 



■In 



\ (l + v / T+4H) 



(7.15) 



(7.16) 



When fe is small — the other limit of interest — the right 
hand side of (7.16) can be expanded, and we find 



\(iu\ 



4- The case of small k 



(7.17) 



Continuing in our investigation of the small- fc regime, 
we note that when k <C 1 we can ignore the difference 
bet ween sinh(fcn ) and sinh(fe(no+l)) and sinh(fc(n +2)) 
in (7.131, and we can also replace sinhfc by k. Then, the 



condition on k is 

k = (l - 4) sinh(/cn ) e - fc (7.18) 
This limit can also be applied to the calculation of the 



properties of bound states. For example, Eq. (7.9 1 re- 
duces to 



tan 9 



-{z/z' - l)sin 2 (fcn )/fc 
1 + {z/z' — 1) sin(fcno) cos(fcno)/fe 



(7.19) 



As will be verified in Section VIII Eqs. ( 7.17 )-( 7. 19 1 



are consistent with the equations at which one arrives in 
the case of a delta function potential in the corresponding 
Schrodinger equation. This is especially the case if we 
expand the factor (1 — z/z') to first order in (3u, which 
corresponds to the magnitude of the Dirac delta function 
potential. 



VIII. DELTA-FUNCTION POTENTIAL WELL 
IN THE VICINITY OF A SURFACE 

Henceforth, we will assume that the functions, 4> n {z) 
can be rewritten in the form <f){x, z) where x is now a con- 
tinuous variable, and that those function can be deter- 
mined by solving the appropriate version of a Schrodinger 
equation. We then assume an attractive delta function 
potential in the vicinity of the bounding surface at x = 0, 
and we rederive the equations satisfied by the bound state 
and the phase shift in extended states. We start with the 
standard "matching condition" 



4>'{x) 
4>{x) 



4>{x) 



V 



(8.1) 



x— x\+e 



where V is the strength of the attractive potential at 
the point x\, which we use as a simpler substitute for 
the combination {z/z' — 1). Note that the dependence of 
4>{x, z) on z has been suppressed. This practice will be 
followed throughout this section. 



A. The bound state 

The unnormalized bound state is given by 
cj) b {x) - 



sinh kx 
sinh nx±e 



c x>x\ ( 8-2 ) 



Applying the matching condition (8.1l to the solution 



2 ) , we end up with the equation for the quantity k 

(8.3) 



coth KXx + 1 



V 

K 



One readily establishes, by looking at small-re limits, that 
there is now a minimum value of V required to sustain 
a bound state. In order for this to be possible, we must 
have V > l/x\. An alternate, but equivalent, version of 



3) is 



re - Ve KXl sinh kxi = 



•4) 



This relation is to be compared with(7.18l 



The normalization of the bound state is obtained by 
taking the integral 



4>b{x) 2 dx 



1 

2^ 



[sinh reccie" 



KXl I 



(8.5) 
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B. Extended states 



Here, the states are of the form 



, sin , fca: sinfei + 9(k)) x < xi , n „, 

(f)(x) = { s . m ^fi \. y " (8.6) 



sin(kx + 9(k)) 



x > Xi 



Note that the form of the extended eigenfunctions is con- 
sistent with the normalization considerations laid out in 
Appendix [Cj The equation satisfied by 9(k) is 

fccotfcxi - fccot(fcxi + 9{k)) = V (8.7) 



After manipulations like those in Section [VII[ we end up 
with the equation for 9{k) 



¥ sin 2 kx\ 



tan0(fc) = 



1 — X- cos kx\ sin fcxi 



This relation effectively replicates (7.191 



C. The calculation of the generating function 

We start by noting that the following holds in the vicin- 
ity of the critical point. 

1 1 



4 cos 2 k 



k 2 



(8.9) 



Given the values that k takes in a system with constant 
potential and large system size, L, (see ( 4.14[ )), the poles 
in the density as a function of fugacity, z, will be closely- 
spaced on the negative z-axis as indicated in Fig. [9] 

Im(z) 



Re (z) 



FIG. 9: The locations of the poles in the density as a function 
of the fugacity, z. The location of z c — —1/4, is indicated. 



We will start with the generating function for the gas 
of rods deep in the interior of a constant potential system. 
We use the Mittag-Leffler theorem [30] to reconstruct the 
density from its poles. In the case at hand, we can ignore 
the spatial structure of the modes, so that the residues 



can be assumed to be constant. The relevant expression 
is 



1 



1 



z 

i 



Zl 



(8.10) 



where quantities zi are the poles and 1/L encapsulates 
the normalization of the modes. Given (8.9 1, this sum 
reduces to 



1 



27r|z r 



dk 



1 



, (z+\z c \)/\z c \ + k2 ^W\ 2 VW^ Z 

(8.11) 

As previously, z c = —1/4. Obtaining the coe fficien t of 



in the expansion of the right hand side of (8.11 1 is a 



pretty straightforward exercise. I will utilize a method 
that is, at least initially, a bit more complicated. We start 
with the standard contour integral-based expression for 
the coefficient of z N in the expansion of a function of z. 
As applied to the function at issue here, it is 



1 



1 



'2m 



1 



1 



dz 



.12) 



where the closed contour encircles the origin. We distort 
the contour so that it wraps around the branch cut that 
starts at z — z c and extends to — oo. If we write z = 
— \z c \ — ( then the integral will look like this, to within 
overall multiplicative constants. 



1 



1 



(-\z c \-C) N+1 VC 



(-1) 



N+l 



1 



1 



/o (\zc\+() N+1 VC 
Now, we introduce the following identity 

1 1 r°° 

1 t N e- t( - z ^dt 



d( (8.13) 



(Iz.l + C)^ 1 
The double integral we now have to perform is 



1 



(-Dm' * 



d(t 



N ' 



-t(W+C) 



'o Jo vC 
We will perform the integral over £ first. It is 



.14) 



.15) 



d( = 2 



dy 



o-ty' 



dy 



(8.16) 



The remaining integral is 



(-1) 



N+l ' 



t N-l/2 er \z c \t dt 



(-1) 



N+ m J ^[{N-l/2)\nt~t\z c \]dt 



.17) 
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When JV 3> 1 , we can evaluate the integral by looking for Substituting this back into the integrand we find for the 
its maximum. The extremum equation that determines result of the integral 
the maximizing value of t is 

^ -1^1=0 (8.18) 
The solution to this equation is 

t=^fi (8-19) 



(-1) 



JV+l 



T M exp 



-1) 



JV+l 



vJV+1 



,A--1/2.1u(^i^)-(.Y -1/21 



"tt— exp [(JV - 1/2) In \z c \ + (JV - 1/2) In JV - 1/2 - (JV - 1/2)] 



JV! 



^-(N-i/2) exp [( N _ 1/2) hi JV-JV] 



-l) JV+1 V7F|^ c r (JV " 1/2) exp [(JV - 1/2) In JV - N - (NUN - N)\ 



^\z c \- {N - 1/2) exp [-(1/2) InJV] 



-l) N+1 V^\z c \- {N - m N-^ 2 



(8.20) 



In the fourth line of (8.201, Stirling's formula for N\ was 
used. 

The principal result, that the coefficient of goes as 



-A' 



JV -1 / 2 , could have been derived considerably more 
easily. However, the method can be generalized. For ex- 
ample, consider the case of the generating function near 
a boundary. Here, we have for the generating function 



sin kx 
\z c \ + z + k 



-dk 



1 — cos 2fcx 



o 

poo 



\z c \ + z + k 2 
1 — cos 2kx 



2tt 7^00 \z c \ + z + k 2 

1 



dk 

dk 



-2y/\z c \+z x 



2^\z c \ + z 



.21) 



We know how to extract the coefficient of z in part of 



the expression on the last line of (8.211. For the addi 



tional part, the double integral corresponding to 
is 



(-1) 



JV+l 



1 

JvT 



dt 



dC t 



JV 



-*(W+C) 



COs(2-\/C x) 



(8.22) 

Again, we perform the integral over £ first, changing in- 
tegration variables as in (8.16). We end up with the 
integral 



-ty 2 



cos(2y x)dy 



,-ty +2iyx 



dy 



.23) 



The final integration to perform is 

{~l) N+1 J o °° exp [(JV - 1/2) In t t\z c \ - x 2 ft] 



The new extremum equation is 
JV - 1/2 



z c + 







(8.24) 



25) 



t ' t 2 

The analysis can be short-circuited if we take into ac- 
count the following facts: 

1. The correction to the solution of the equation due 
to the last term will be small. 

2. The effect of the correction on the first two terms 



in the exponent in (8.24) will also be very small, as 



t has already been adjusted so that those terms are 
at an extremum. 



15| This all means that the re sult o f the integration in (8.241 



will to be the same as in (8.171, except that there is the 
additional term 



t 



N -1/2 



.26) 



Combining this with the term we have already evaluated 
we have for the coefficient of z N in the case of the hard- 
rod gas near an end-wall 

(-1)^+1 V^cl-^-^JV" 1 / 2 (l - e -\^/iN-l/2)^ 

(8.27) 

The difference between JV and JV + 1/2 can be neglected 
in the denominator in the exponent in (8.27). 
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D. Density in a finite interval 

We can also utilize the Mittag-Leffler method to recon- 
struct the density in the case of a finite interval. Here, 
the reconstructed density is, to within an overall multi- 
plicative factor 



p{n) 



\t(\ 



L ^ V 4 

m— 1 
AN 00 



4i 2 



sin 2 (nmn / L) 



E 



-NHl + ^^/^) sin 2 {nmn/L) 



m— 1 

<N 00 



L 



^2 e~ w m N/L sin 2 (nmn/L) 



(8.28) 



p(n)4 




10000 



n 



m—l 



FIG. 11: The density as given by the last line of (8.28 1, for a 
gas of hard rods, normalized as in Fig. [10] In this case the 
values of N are 50, 5,000, 500,000, 5,000,000 and 50,000,000. 
Again, the highest curves correspond to the smallest values 
of N. 



where the explicit value of z c is used. There are two 
different limits to consider, based on the ratio N/L 2 . If 
N ^S> L 2 , then the sum is dominated by the first term, 
and the density is 



p{n) = e -*~ N,L sm 2 (irn/L) 



.29) 



On the other hand, if N <C L 2 , then the sum is as given 
by (p7l. 



Figure 10 shows how the density as given by (8.28) 



behaves as a function of n when L = 10, 000 for various 
values of N. The function is multiplied by 4r N \fN so 
that the various curves tend to the same value in the 
interior of the interval when N is small enough. 



p(n)4 




n 
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FIG. 10: The density, as given by the last line of ( |8.28[ l, as 
a function of n, multiplied by A~ N y/N, with L = 10,000, for 
the following value of N: 50, 500, 5,000, 50,000. The highest 
curves correspond to the smallest values of N . The plot is 
restricted to the region adjacent to the boundary. 



Figure 11 displays the density profile over the entire 
interval for a different set of values of N . The domi- 
nance of the single Schrodinger equation eigenfunction, 
sin(7r(n+ !)/£), is evident in N = 50,000,000 curve. 



E. The attractive potential 

We can make use of the results above to perform the 
integral needed to reconstruct the density as a function of 
the fugacity. An important precursor to the calculation of 
the generating function is the reconstruction of the sum 
J2i 0;( x ) 2 - The details of this calculation are contained 
in Appendix [D] Making use of the results of the results 
of that appendix, we find that the modification of the 
generating function for the density due to the presence 
of the delta function potential near the boundary is 



sin 2 kx\ Re 



k — V sin kxie lkxi 



dk 

(z+\z c \) + k 2 

sin 2 KX\e~ 2K,{ - x ~ x ^ 



h. I sinh ' 



2^ lo.ia.i.a1 K,X\G 

2/ 



sinh 2 ( v /z+ \z c \xi) 
V z +\ z c\ 

e -2^/z+|z c | xy 



y/^+\z^\-V S il^{y/I+\z^X 1 )e-y^+^^ 

(8.30) 



In the above equation, the quantity k satisfies (8.3) or, 
equivalently, (8.4 1. This contribution to the generating 
function is, recall, in addition to the contribution that 
one derives in the absence of the attractive potential. The 



results in (8.301 are relevant to the case x > x±. Now, 



the extraction of the actual density at fixed monomer 
number, TV, 



frc 



30 1 entails the kind of contour inte- 



gral described in Section [TTJ The result of that integration 
depends on the value of V. If the potential strength is 
sufficiently great that there is a bound state, then one can 
show straightforwardly that there is a pole in the last line 
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of (|8.30| with a residue that yields the normalized bound 

-JV 



state with a prefactor going as z 4 , where 



\J\z c \ + z„ = K 



.31) 



The quantity k, again, satisfies the equivalent equations 
3) or, equivalently (8.4), for the bound state. Note 

\z c \, in that 



that the absolute value of z* is less than 

= K 2 — | z c | . We are assuming a k that is not too large 
If V does not exceed the threshold for a bound state 



then things are a bit different. There is no pole in (8.30 1, 
but rather a branch cut. A detailed calculation, under 
the assumption x ~ vTV and x\ <C VN, yields 



p{x) 
oc \z c \ 

1 



-N N -l/2 
_ e -\z c \x 2 /N 



2V\z c \xix/VN 
1 - Vxi 



.32) 



Note the denominator 
the left hand side of 
ground state. 



in the last term in brackets on 
.321, a signature of an impending 



is multiplied by z^ , where is the value of z correspond- 
ing to the pole associated with the lowest energy bound 
state. In the case of Fig. [12] the attractive potentials are 
sufficiently strong to ensure two bound states. Figure fl3] 
displays the same set of modified densities, this time over 
the entire interval. 




10000 



FIG. 13: The same set of densities graphed in Fig. |12[ this 
time over the entire interval. 



F. Attractive potential at the boundaries of a finite 
interval 



Some numerical results serve to illustrate the effects 
of attractive potentials on the density of branched poly- 
mers confined to a finite interval. For example, Fig. [12] 
shows what the density looks like for various values of the 
number of monomers, TV, when there is an attractive po- 
tential a distance n — 40 from the boundaries of a region 
with extension L — 10, 000. Here, we take TV — 50, 500, 
5,000 and 50,000. The figure graphs the density close in 




n 



100 



200 



300 



400 



FIG. 12: The density, p(n), for the case of attractive po- 
tentials near the bounding surfaces of of a region of length 
L = 10, 000. The curves correspond to TV =50, 500, 5,000 
and 50,000. The heights of the density curves decrease with 
increasing TV. 



to one of the bounding surfaces. In each case, the density 



Finally, we consider the case of an attractive potential 
that is not quite strong enough to generate bound states. 
Here, one might expect to see the influence of a "precur- 
sor effect." Figure [Ti] graphs the density, normalized as 
in Fig. |10| The dependence on number of particles, TV, 
which is not as simple and monotonic as in the previous 
cases, is indicated in the figure. The dependence of p(n) 



50000 




100 200 300 400 



n 



FIG. 14: The density, normalized as in Fig. |10| in the vicinity 
of one of the boundaries of a system with length, L — 10, 000, 
for TV equal to 50, 500, 5,000 and 50,000. The attractive 
potential, which is not quite strong enough to induce a bound 
state, is at a distance 40 from the boundary. 



on TV in this case is to be contrasted with the graphs of 
p(n) in Fig. 



10 



where at TV = 5, 000, the normalized 
density is greatly suppressed with respect to the density 
at lower values of TV. 
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IX. THE PRESSURE OF ROOTED TREES 

Recall that the generating function for the number of 
trees with roots at the position x is derivable from the 
density of a gas of hard-core particles. As we have seen, 
the form of this density is 



E 



z k + z 



(9.1) 



The functions ipk(x) arc solutions to a Schrodingcr-likc 
equation that are, furthermore, normalized. If we sum 
over all possible locations of the roots, we end up with a 
generating function directly derivable from the sum for 
all rooted trees in the interval. This generating function 
is 



E 



i 



Zk + z 



(9.2) 



We are interested in the coefficient of z in this sum, 
which leads us to the following analogue of the partition 
function for rooted trees in the interval: 



E- 



-JV 



(9.3) 



There are additional combinatorial factors, but, as we 
will be looking at large N and taking a log, they turn 
out to be unimportant. 

From the expression above, we obtain the following 
result for the effective free energy: 



k B Tln(^2\z k \- N ^j 



(9.4) 



If N is large enough that the z k closest to the origin 
dominates, then the free energy reduces to 



Nk B Tln Zl 



(9.5) 



where Z\ is the location of the singularity that lies closest 
to the origin. 



A. Pressure in the absence of a potential energy: 
scaling formulas 

In the case of an interval with no potential, we can 
write 



:/, = - - I t J 



(9.6) 



where z c — w = 1/4. Then, the free energy has the form 

-N~ 



fegTln 



Eh + 



w 



NksTln z c — ksThi 




-N 



(9.7) 



The last line of (9.7) follows if L is sufficiently large. We 



assume this to be the case and proceed. The pressure 
is the negative of the derivative of the free energy with 
respect to L. Making use of (9.7) we find 



P = 2k B T"~ ^ k L 



_ N ^ {kw/L) 2 



(9.8) 



Keeping the full expression for the free energy, we ob- 
tain a general result for the force between the two walls. 
Figure [15] illustrates the dependence on L of the force 
divided by k B T, with w — z c — 1 and N = 100. One can 

P/k B T 
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FIG. 15: The pressure as a function of L when N = 100. 
Here, w — z c — 1. 



easily show that at sufficiently large L the pressure will 
go as 1/L. 

In fact, a straightforward analysis tells us that the pres- 
sure will, in this case have the general form 



P(L, N)/k B T = V(L/VN)/L 



(9.9) 



Figure |16| is a graph of the function on the right hand 
side of (9.9 1. The large L behavior of the quantity V is 



evident from the figure, in that it approaches a limiting 
value as its argument goes to infinity. 

There is also the question of the behavior of the pres- 
sure at intermediate values of L\J~N . In the regime 
in which one singularity in z dominates, one can show 
that the pressure goes as N/L 3 . This is consistent with 
V cx N/L 2 . Figure fTT] is a plot of (L/VN) 2 x V. 
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FIG. 16: The quantity LP(L,N)/k B T = V (L/y/N). 
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FIG. 17: The quantity (L.y/NfV, plotted against L/y/N. 



Note that this combination is effectively a constant 
between L/y/N = and L/y/N = 3. From this we can 
infer two different regimes for the pressure. In the first, 
when L ^> y/N, the force goes as 1/L, independent of 
N. In the second in which L/y/N ~ 1, the force goes as 
N/L 3 . 



X. THE INFLUENCE OF ATTRACTIVE 
POTENTIALS 

Once again, we assume that there are attractive po- 
tentials near the two bounding surfaces, as indicated in 



Fig. 18 There will be two types of cigcnfunction in this 
configuration: even and odd parity about the center. In 
the region < y < L/2, the even functions have the form 

/ ( \— f sm ^2/ U< x 
V^y* ~ \ cos(fc(L/2 -y)) x<y< L/2 ( > 

The odd parity functions are of the form 

^"^ ~ \ sin(fc(L/2 -y)) x<y< L/2 ^ U '^ 

The equation for the eigenvalues, expressed through the 
quantity k, is 



ip'(x-) 4>'(x+) 



il>(x~ 



V 



(10.3) 



FIG. 18: The potential configuration with which we will work. 
There are two attractive delta function potentials, each a dis- 
tance x from the bounding walls of a region with an extent 



where V is the strength of the delta function potential, 
x~ is just to the left of x and x + is just to the right. 
Making use of the two forms in (10.1 1 and (10.21, we 



have for the equations satisfied by k in the case of the 
even and odd cigcnfunctions, respectively 

kcot(kx) + ktan(k(x J J = V (10.4) 



V (10.5) 



kcot(kx) — A: cot fc I x 



For sufficiently large values of V, there may also be so- 
lutions at imaginary k = in. The equations that are sat- 
isfied in the even and odd parity cases are, respectively, 



KCoth(K2;) — fttanh ^x — 
Kcoth(ra) — Kcoth I k I x — 



= V (10.6) 
= V (10.7) 



The threshold values of V for which there arc solutions 



to (10.61 and (10.7) are, respectively 



^cvcn 
V dd 



1 

X 

1 



1 



L/2-x 



(10.8) 
(10.9) 



Given these solutions, we are able construct the ex- 
pression for the force generated by the rooted branched 
polymers. As previously, we start with the result (9.4 1 



for the effective free energy. The general result for the 
force associated with this free energy, corresponding to 
(pOl is 



P 

ksf 



d_ 

OL 



\-N 



-NY, l \z l \- N - 1 d\z l \/dL 



-JV 



(10.10) 
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To assess the derivative of the singularities with respect 
to L, we make use of the general relationship 



P/k R T 



i 



Then, 



d\zi\ 
dL 



4 cos 2 ki 



sin ki dki 
2 cos 3 ki dL 



(10.11) 



(10.12) 



There will be corresponding equations for the up to two 
bound states, in which an analytic continuation has been 
performed from the variable k to the variable in. 



Equations ( 10.10 1— (110.12 1 allow us to calculate results 



for the force exerted by the branched polymer between 
two walls. Figure [19] summarizes results for the pres- 
sure for various values of the attractive potential. In 
the case of interest here, there is no bound state un- 
less V > 10. Note that for sufficiently small separa- 



P/k B T 




-0.002 



FIG. 19: The pressure as a function of distance, L, between 
the walls for x — 0.1 and the following values of the attractive 
potential strength, V: 9.9, 10.01, 10.02, 10.05. The smaller 
V the higher the pressure. 



tions and sufficiently attractive potentials the pressure 
becomes negative, corresponding to an attraction be- 
tween the walls. This follows from the fact that the 
bound state free energy decreases as the attractive wells 
approach each other, in analogy to the simplest version 
of the chemical bond [31] . 

Figures |20] and [21] are two interesting and contrasting 
plots. In both figures, the attractive potential is a dis- 
tance 0.1 from the edges of the system. The threshold for 
a bound state isV— 10. In the first plot, V = 10.00025, 
which is just sufficiently strong that there is a bound 
state In the second plot, V = 9.9997, and the attractive 
delta function does not quite suffice to produce such a 
solution to the Schrodinger-like equation. The plots are 
for a range of values of TV, as indicated in the caption. In 
the case of the first plot, Fig. [20] the larger the value of 
TV, the lower, or more negative, the pressure. In the case 
of the second plot, Fig. [21] the greater TV, the higher the 
pressure. 

What this tells us is a very small change in the attrac- 
tive potential suffices to give rise to a considerable change 
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FIG. 20: The pressure divided by k B T for V = 10.00025 and 
the following values of TV: 300,000, 100,000, 30,000, 10,000, 
3,000. The larger TV the more negative (or less positive) the 
pressure. 
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FIG. 21: The pressure divided by k B T for V = 9.9997 and the 
following values of TV: 300,000, 100,000, 30,000, 10,000, 3,000. 
Here, the larger TV the greater (more positive) the pressure. 



in the force between the walls. It also tells us that an 
attractive potential can be "mediated" by the branched 
polymer in such a way as to facilitate the assembly of a 
capsid. 

Finally, note in the case of Fig. [20 that the force is con- 
cave downward when L is small enough. This leads to a 
mechanical instability in that range, in which the walls 
will continue to collapse, assuming a constant counter- 
vailing force. This means that the walls will be pulled 
together until other mechanisms, such as excluded vol- 
ume effects, intervene. 
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APPENDIX A: NORMALIZATION OF 
EIGENFUNCTIONS 

We start with the two equations 



dx 2 
d 2 fo{x) 
dx 2 



(Al) 



+ V{x)cj)2{x) = (E + AE)<h(x) (A2) 



We will assume that 4>i{x) has zeros at x = and x — 
xq. As for 02 (x), it is also zero at x = 0, and it has 
another zero close to xq. In fact, we assume that as 



AE -> 0, <j> 2 {x) -> 4>i{x). Subtracting (|AT) from ( |A2| ) 
and integrating from x — to x — Xq, we find 



f 

Jo 



d^fnix) ,d 2 cj} 2 (x) 
^ >~ dx 2 ^ ' 



dx 2 



dx 



dx 



.d(j>x{x) ,d(j> 2 (x) 
Mx) ^x-~ Mx) ^x~ 



dx 



dfaix) 

$2 \Xj — 01 (X) 

dx dx 



AE I <j) 1 {x)<j)2{x)dx 



(A3) 



As for the next to last line in ( A3 ) , because of bound 



ary conditions all contributions are equal to zero except 
4>2(xo)(d(f>i(x) / dx)\ x=Xo . Under the assumption that AE 
is small, the last line reduces to AE <fri(x) 2 dx. That 
is, we now have 



02(x o ) 



d(j>i{x) 



dx 



AE I ^{xfdx (A4) 
la 



Given that the zero of <p2 (x) is close to Xq, we can write 



4>2{xq) 



d(j>i(x) 



dx 



dxQ 
IE 



AE 



(A5) 



This can be established graphically. Combining the re- 
sults, we have 



dcj)i(x) 



dx 




Now, let us look at the denominator 0„0„_i in (4.371 
In the vicinity of the zero of </>„ , we can write 



where 



(n - x (z)) 



( \ , dXn A 

dz 



(A7) 



(A8) 



On the other hand 



This means that we have 



b n -i = 4>{— 1 — Az dxo/dz) 

d(j>n 

dn 



—Az dxo/dz)(f)(—l) 

2 



d(j) n \ dx 



dn 



dz 



Az 



(A9) 



(A10) 



If we divide by this, we have a contribution to the residue 
of the pole going exactly like l/(/ 4>n dn), to within a 
multiplicative constant. In fact, a careful analysis of the 
relationship between the variables z in our system and 
the energy, E, in the Schrodinger equation leads to the 
conclusion that the multiplicative constant is precisely 
1/4. This means that, to within the multiplicative factor 
of 1/4, the residue effectively normalizes the contribu- 
tions to the density. 



APPENDIX B: NORMALIZATION IN THE 
DISCRETE SYSTEM 

We start with the following version of the equation for 
the quantity 4>n 

e -/3«n+l/2^ n+l(fc) + e -/3«n/2^_ l(fe) = 2cQsfc 0n(fc) 

We consider two versions of this equation, one just like 



(Bl ) and one with the parameter k slightly different. De- 



noting that new value of k as fc', we have, multiplying 
(Bl ) by (j> n {k') and the corresponding equation for k = k' 



4>i(x) dx (A6) by 4> n (k), subtracting the results and summing over n, 
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£ Mk') (e-^/ 2 <p n+1 (k) +e- /3 "«/Vn-i(fc)) 

i=n b 

- £ 0„(fc) (e"^ 1 / 2 ^ +1 (fc') + e-' 3tl " /2 0„-i(A ; ')) 

n—nj, 

- £ (0„(fc)e-' 3 ""+ 1 / 2 ^ +1 (fc')-^(fc')e- /3 "" /2 0„-i(fc)) 

n— nfa 

-0nMe~^ +l/2 ^ + i(k') + K(fc')e^"^ /2 ^-i(fc) 

7T. e 

= 2(cosfc-cosfc') ^ (j>n(k')<f>„{k) (B2) 

n — n 5 

i 



We will take the solutions of the equations in the equation 
above to satisfy the boundary conditions in Section |IV| 
Furthermore, we take the end point of the summation to 
be n& = 1, n e = L — 1. Furthermore, we will assume that 
the potential energy is zero at and near the boundaries. 
Then, (f) rLb {k) = 4> nb -i{k) and similarly for k' . Addition- 
ally, we will assume that the value k is consistent with the 
boundary condition <f) nc =L{k) = 0. Then, (B2| reduces 
to 

L-l 

- 0£_i(*r)0 z (ifc / ) = 2(cos k - cos k') ^ </>„(fc)0 n (fc') 



A. 



(B3) 

From here on, the analysis follows that in Appendix 



Given that, in our versi on o f the Schrodingcr equation, 
E = k 2 , we can rewrite (C3) as 



dxa 1 dxo 
~d~E ~ 2k~dk 

2k 2 



From ( A6 ) and the above, we have 
V dx J X=XQ 



x + 0'{k) 
2k 2 



(C4) 



(C5) 



We now return to (CI I. If the integer I increments by 
one, then k will change as follows 



Ak(x + e'(k)) = 7T 



(C6) 



APPENDIX C: NOTE ON THE 
NORMALIZATION OF EXTENDED 
EIGENSTATES 



This tells us that 



(C7) 



As an essential step in the calculation of generat- 
ing functions, we establish the proper way to normal- 
ize the eigenfunctions we deal with in the case of very 
large systems. In particular, we are interested in the 
case of an eigenfunction in a long interval that goes as 
<j>{x) oc sin(kx + 8(k)) towards the left end of the interval. 
The boundary condition will be 4>{xq) = 0, which tells us 
that 



kix Q + 9{ki) = In 



Recall ( A6 ) . For a given value of ki , we have 
Ax(x + 9'{h)) + x Ak = 



This tells us that 



dx 
~dk 



x 



9'(k) 



(CI) 
(C2) 

(C3) 



and, we then have the following result for the eigenfunc- 
tion sum leading to the density: 



E 



f*° 4>i(x) 2 dx 

E^ 



2kf 



Akixo + e'ih)) 



(d<f>(x)/dx)U X0 fa + O'ih)) 
k 2 



6{x) dk 

*J {dcp{x)/dxf x=x J K 



(C8) 



If we choose 4>(x) precisely equal to sin(fcx + 0(k)) at x 
near xq, then the last line in (C8l reduces to 



- / <j)(x) 2 dk 

7T 



(C9) 
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APPENDIX D: RECONSTRUCTION OF THE 

sum E,<Mz) 2 - 

We focus on the case x < X\ and begin with the ex- 
tended states, and we will seek the difference between 
this sum and the sum in the absence of the attractive 
potential. That is, our task is to find the value of the 
sum. 



Esin kx r . 2 . 
— 2 sin (kx\ 
sin kx\ 



+ 9(k)) -sin 2 (fcri)] (Dl) 



After expanding sin(foci + 9) and carrying out some sim- 
ple algebra and trigonometry, the above expression re- 
duces to 



E- 



l kx 



sin kx\ 



D 2ikx! / ( 



,2i6(k) 



1 



Given the trigonometric identity 



The summand becomes 
1 sin 2 kx 

Re 



2 sin kx± 



sin kx Re 



i tan 9 + 1 
1 — i tan 9 



_2ik Xl 2 »tan6>(fc) 
1 - itan0(fc) 

^ e 2ikx! y 

k — V sinkxie lkxi 



(D2) 



(D3) 



(D4) 



where (8.8) has been used. The integral to be performed 



IS 



1 



sin 2 kx Re 



,,2ik:i 



V 



dk (D5) 



k — V sinkx\e lkxi 
where the removal of a factor of two in front of the in- 



tegral (see (C9 1) is compensated for by the fact that 
the range of integration has been extended by a factor 
of two. In fact, the result of the integration will be 
real as a matter of course, so we can remove the "Re" 
function from the expression. The sole contribution to 
the integration is from a pole on the upper imaginary 
axis, at a value of k = in, such that the denominator 
k — V sin kxie lkxi — > i(k — V sinh nxie~ KXl ) = 0. This 
equation for the pole is the same as the requirement 
for the bound state. The residue at that pole is 



4) 



iV sinh 2 Kxe 2kXi 



— x 2iri x — ; 
7T - V smkxie lkx i)) k _ 



(D6) 



Focusing on the denominator in (D6), we have 
d 



dk 



(k-Vsmkxie llcXl ) 

(l-Vx ie ^) k=iK 
1 - Vx ie - 2KX1 



= l - 



sinh kx 



sinh kx\ 

KX\ 

(sinh nx\e KXl 



KXl) 



(D7) 



The third line of (D7l follows from (8.4 1. Inserting this 



result into (D6l we end up with the result for the inte- 
gration 



2Ve KXl sinh nx\ . 2 

sinh kx 



sinh 



KXl 



2 k 



sinh 



sinh kx 



KX\e K 



KX\ 



(D8) 



Again, we have utilized (8.4 1 to obtain the right hand 
side of the equation above. Referring to Section [VIII A| 
we see that this exactly cancels the contribution of the 
bound state to the left of the delta function potential. 
We have recovered the standard result for a complete 
orthonormal set of cigcnfunctions, and that is that the 
sum is independent of the set chosen. 

The corresponding integration when x > x\ is 



sin kx\ Re 



„2ik 



xy 



k — V sin kxie lkxi 



dk (D9) 



The same sort of contour integration yields a result that 
precisely cancels the contribution of the bound state in 
that regime. 



APPENDIX E: DIMENSIONAL REDUCTION 
FOR A SPHERICAL GEOMETRY 

A prime motivation for the work reported here is the 
packing of complex RNA into a viral capsid. The closest 
approximation to the geometry of the packing environ- 
ment entails spherical symmetry. It is natural to ask 
whether dimensional reduction will prove useful in this 
case. In this Appendix, we explore the consequences of 
the results of Brydges and Imbrie [17, 18 when spherical 
symmetry holds. 

The two physical systems related by dimensional re- 
duction are a hard-core classical gas in D dimensions 
and a solution of rooted branched polymers in d = D + 2 
dimensions. The particles have a hard core repulsive in- 
teraction with diameter a, and they are confined in a 
container with radius R. The equation relating these 
two systems is: 



Ghg(z) = G BP {- 



2-na 2-kR 



(El) 



where Ghg on the left side of the final equation is the 
grand partition function of the gas, z being the fugacity 
of the particles. The quantity Gbp(z,w), on the right 
side of the equation, is the grand partition function of 
the solution of annealed branched polymers. The first 
argument, z, is the fugacity of the monomers and the 
second argument, w, is the fugacity of the roots. The 
roots arc constrained to be on the surface of the container 



(see figure (22)) 



We provide the proof of the equation (Ell in the fol- 



lowing sections. The reasoning is parallel to the one sug- 
gested by Cardy [13] • We start from the partition func- 
tion of hardcore classical gas in. Then we continue with 
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If we define: 




FIG. 22: Solution of branched polymers. Roots are depicted 
by solid black circles. 



defining the partition function of a super gas. We show 
those two are equal using cluster expansion and the prop- 
erties of gaussian integration. Finally we put the parti- 
tion function of a super gas into the form of the partition 
function of a solution of branched polymers. This last 
step is done using the Taylor expansion of functions of 
Grassman variables. 



and 



9rj = 9rl = exp(L7"Cy)) ~ 1 



hi = h r 2 = cxp(V(r^)) 



then it is well known that the grand partition function 
can also be represented in the form of cluster expansion 
[32] - 



(E4) 



Ghg(z) = exp(J2 b k zk ) 

k 

k 

b " = ^E/ U^nhiJlgij (E5) 

(E6) 



eg 



The label CQ stands for 'Connected Graphs'. is the 
sum of all connected graphs (clusters) with k particles 
divided by k\. in the next section we define the partition 
function of a super symmetric gas. 



2. super-symmetric gas 



1. Classical Hardcore Gas 

Consider the classical gas in a D dimensional space. 
We assume the interaction potential energy of two gas 
particles, U(rfj), depends only on the distance between 
those two particles: rfj. The labels, i and j, indicate 
the gas particles. We also include an external potential, 
V(rf). This potential depends on the radial distance 
of the gas particle from the origin: rf. The partition 
functions of the hardcore gas is 

Ghg(z) = Y, Z hg{N)z n (E2) 

Zhg(N) = ± jlK'v >: ' " r ' ' 

i 

(E3) 



Consider the supersymmetric gas in a super-space of 
d = D + 2 real plus 2 Grassman coordinates. As in 
Cardy's exposition [19] the distance between two parti- 
cles and also a particle from the origin is defined as: 

R% = (r<-rO 3 + (0<-0i)(0i-0i) (E7) 
R? = tj + BiOi (E8) 

The quantities and j-j are the usual real-valued dis- 
tances in d dimensional real space, while and 9i are the 
Grassman coordinates of the particles. As in the case of 
the classical hardcore gas, it is natural to define the par- 
tition function of a supersymmetric gas as ('SSG' stands 
for 'Super Symmetric Gas'): 
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Gssg(z) = Y. Zh g{N) 
1 

m 



,N 



^SSG 



(N) 



■T. ij U(R%)-Y. i V(R 2 i ) 



(E9) 
(E10) 



r 



Here, the cluster expansion works as well. The only modi- 
fication is that wherever we have or rf in the classical 
case, we use their counterparts for the supersymmetric 
and R 2 . In the next section we demon- 



case, 



R ij 



i.e. 

strate that, given a supersymmetric gas in super-space 
of d — D + 2 real and two Grassman coordinates and a 
hardcore classical gas in a D dimensional space with the 
same inter-particle and external potentials, we have 



G 



SSG 



(z) = G HG (z) 



(Ell) 



The two Grassman coordinates have canceled the effect 
of two real coordinates. It is evident that if we start with 
a classical hardcore gas in D = d — 2 dimensions, we 
arrive at the same expression as the above. 



This completes the proof for the equation (Ell ) 



In the next section, we expand the partition function 
of a supersymmetric gas. However, instead of cluster ex- 
pansion, we use a Taylor expansion. It turns out that the 
super-symmetric gas is related to the generating function 
of branched polymers. 



3. Super Symmetric Gas versus Classical Gas 

As in [T5], we define the functions p(/j,) and q(v) as 
follows: 

r°° 2 

g(rl) = / <i/',,/": /',,» " ' 



K r i) = / duiq(ui)e~ 
Jo 



4. Super Symmetric Gas versus Branched 
Polymers 

We begin this section by defining four new functions: 

P(4) = e~ u ^ (E12) 
S(r 2 ) = e~ v ^ (E13) 



Notice that if we Taylor-expand P(R 2 j) and R(R 2 ) 
Notice that in the above expressions, R% can also be around r l and r ] > because of the Properties of Grass 
used instead of r 2 4 . Consider the contribution of of a 
connected graph of k supersymmetric particles, CQ. We 
use the above relations to transform the functions g and 
h into the functions p and q. Regardless of integrations 
over the parameters jiij and Ui, we are left with 

k 

f J] d d rM exp(-i ]T {r Af n + S^fr) 
J i=i ji 

The quantity A c ^ is a k x k matrix dependent on the con- 
nected graph CQ. The gaussian integration over Grass- 
man numbers evaluates to 



{- l -k) k &et{A c S) 
We have used the following convention in this evaluation: 
d£d(9 66 = A 



On the other hand the integration over the real coordi- 
nates contribute a factor of 

(2n) kd / 2 det(A CG )- d / 2 

If we choose A = — - , the final expression evaluates to 

(27r) fe(d-2)/2 det(A Cff ) ( ( i-2)/2 



man numbers, only the first two terms survive 

P(R%) = Ptfj+OydijQtfj) (E14) 
S(R 2 r ) = P(r 2 ) + 6AT(r 2 ) (E15) 

In the above equations, Q and T are the first derivatives 
of the functions P and S respectively. Using the above 
equations, we can expand the grand partition function 
of a supersymmetric gas. We can construct a graphical 
expression for each term. We indicate the term Q(r^) 
by a pair of connected particles (monomers) at ad Tj. 
The term T(rf) is represented by a solid black circle at 
point Ti representing a root monomer. For each non-root 
monomer, represented by a hollow circle, located at 
we multiply a factor of S{rf). For any pair of monomers, 
located at Vj and rj, which are not connected we multiply 
a factor of P(r?-). For each monomer we also have a 
factor of z. The supersymmetric grand partition function 
is the sum of all these terms (graphs). Each term has a 
product of its connected graphs. It can be shown that 
any connected graph which has a loop or no root or even 
more than one root is zero (See figure 23 1 . 



Thus, only those graphs which are products of con- 
nected, loopless and single rooted graphs, ie: branched 
polymers, contribute. As the result of Grassman inte- 
gration, for each monomer a factor of — 1/tt is generated. 
In order to go further we need to specify the potentials 
t/(r?.) and V(rf). 
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rfj < a 2 

1 4 > « 2 

' >R 2 







oo r' 



1 



<i? 2 
>i? 2 



Because of this, Q and T take the form of delta functions 

Q(rfj) = Stfj - a 2 ) = ±5{r l3 - a) 

T(r 2 ) = S(r 2 -R 2 ) = -J R S(r i -R) 

With the above choice of functions, any pair of con- 
nected monomers have a fixed separation of a, the di- 
ameter of monomers. Also, any root monomer, because 
the delta function is constrained to stay on the surface 
of the sphere. For any root monomer we have a factor of 
— 1/2R and for all other monomers a factor of I /2a. Re- 
call that there is also a factor of —z/n for any monomer. 
Therefore, z/2ttR appears as the fugacity of the roots 
and —z/2i:a as the fugacity of other monomers. Putting 
everything together we obtain 



Gssg( z ) = Gbp{- 



2na 2nR 



(E16) 



Comparing this with our previous equation (Ell I pro- 
vides the main result: 



FIG. 23: Some of the vanishing terms in the expansion of 
partirion function 



Ghg(z) = G BP {- 



2-na 2-kR 



(E17) 



Discussion 



We choose a hardcore repulsive potential for U, a being 
the dimameter of monomers. Imagine a spherical con- 
tainer for the particles of radius R. The external poten- 
tial is zero inside and infinitely large outside. This results 

in step functions for P and S: U(r 2 ^) ■ 



oo rfj < a 2 
rfj > a 2 







As we see the fugacities in G bp appear with opposite 
signs in the equation (El I. However, the physical region 



of the partition function of Gb p is where both of its fu- 
gacity arguments are positive. Consequently, one cannot 
explore the physical region of the branched polymer so- 
lution using this equation. 
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